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Charmonium from Lattice QCD 
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Charmonium is an attractive system for the apphcation of lattice QCD methods. While the sub-threshold 
spectrum has been considered in some detail in previous works, it is only very recently that excited and higher- 
spin states and further properties such as radiative transitions and two-photon decays have come to be calculated. 
I report on this recent progress with reference to work done at Jefferson Lab. 
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1. Introduction 



Between 3 and 3.7 GeV a number of states exist 
which are beheved to be the bound states of a charm 
quark and an anti-charm quark and whose widths 
are narrow owing to their being below the thresh- 
old to decay to a pair of open charm mesons cou- 
pled with the suppression of annihilation channels at 
this high mass scale. Because the hadronic contribu- 
tions to their widths are so small, radiative transi- 
tions between them constitute considerable branching 
fractions, and the rates of these transitions have been 
measure d with some accu racy by a number of exper- 
2QQ6| ). Additionally the C = + 



iments (Yao et al. 



states can decay to a pair of photons - this process 
when time-reversed can serve as a production mecha- 
nism (two-photon fusion) at e+e~ machines. 

Rates for these radiative processes have been com- 
puted in various varieties of quark-model, and are typ- 
ically fairly successful when one sets parameters us- 
ing the experimental spectrum, however corrections 
beyond approximations like non-relativistic dynamics 
are often uncontrolled in these models. 

In the current century the charmonium picture has 
filled out considerably and new mysteries have arisen 
owing to the high statistics and new production meth- 
ods made possible by CLEO-c and the 5-factories. 
The remaining expected sub-threshold states, r^^, hc^ 
have been observed, as have radiative transitions from 
the ?/^(3770) down to the XcJ- Th e above-threshol d 

2006] ), 



spectrum is rapidly being mapped (Swanson 



with some states living up to the expectations of po- 



tential models (Uehara et al. 



ing as something of a surprise ( Choi et aLJ 2003] ). The 



2006 ) and others com- 



increasingly complete set of exclusive data in e+e" 
looks set to allow determination of the vector spec- 
trum with some confidence. 



In a series of recent works (Dudek et al. 2007 2006 



Dudek and Edwards 2006 ), members of the Jefferson 
Lab lattice group have investigated the possibility of 
computing excited spectral and radiative quantities 
using lattice QCD. These initial studies have been 
carried out on quenched lattices with rather promis- 
ing results. In the sections that follow I will briefiy 
describe the work done. 



2. Excited and higher spin states 

The mass spectrum of a field theory considered in 
Euclidean space-time can be extracted from the time- 
dependence of a two-point correlation function. 



where Oij are operators that have the right quan- 
tum numbers to produce a particular state from the 
vacuum which are constructed from the fundamen- 
tal fields of theory. For example in QCD we might 
try to study the pseudoscalar spectrum by considering 
an operator ipj^ip. The correlator receives contribu- 
tions from all states in the theory with the appropriate 
quantum numbers. 



ya* ya 



exp —rriat. 



(1) 



(In a quantum mechanical bound state model we 
might think of radial excitations being labeled by 
a). In practice extracting anything other than the 
ground state mass from fits to the time-dependence 
of a single correlator is difficult and often unstable. 
This is particularly troublesome in a system like char- 
monium where there are significant approximate de- 
generacies, e.g. the '0(3686) and '0(3770). These de- 
generacy problems are made worse on a cubic lattice 
where states are labeled not by a continuum spin, but 
by an irreducible representation of the cubic group. 
The continuum spin content of these various irreps 
is shown in Table I. This indicates that, for exam- 
ple, componets of a 3 state would appear in the 
same correlators as 1~~ states. Since from poten- 
tial models we expect there to be a 3~~{^Ds) roughly 
degenerate with the 1 (^^i) '0(3770), we anticipate 
that there should be three roughly degenerate excited 
states above the ground state in a T^ correlator. 
Extracting this from a fit to a single correlator is not 
practical. 

Given this one might consider more reliable ways 
to extract the excited state spectrum. A variational 
method utilizing a large basis of operators satisfies 
this need. Its major advantage is that it utilizes the 



Proceedings of the CHARM 2007 Workshop, Ithaca, NY, August 5-8, 2007 



A dx 



Al 


1 


0,4,6,... 


A2 


1 


3,6,7,... 


E 


2 


2,4,5,... 


Ti 


3 


1,3,4,... 


T2 


3 


2,3,4,... 



Table I The table shows the single- valued irreducible rep- 
resentations A of the cubic group O, together with their 
dimensions c/a and continuum spin content J. Additional 
superscripts are employed to denote charge conjugation C 
and parity P. 



li sm T47i am 

ai X V ^l ai X V sm Qi X V ^^, ai x V J^-f, 

p X P s^t^l p X ID "j^ p X ID s:ii 

TT X B sm 

I :rf21 X B sm 



px©,' 



J- 



•^^^ 



orthogonality of states in a space of operators - while 
states might be degenerate and hence hard to sepa- 
rate on the basis of mass, they remain orthogonal and 
hence easier to separate on the basis of their state 
vectors. 



In |Dudek et al.| [2QQ7| , an operator basis was con- 
structed based upon operators that in the continuum 
would have the structure of fermion bilinears with a 
number of symmetric covariant derivatives 



O 



fiup-- 



= ^{x)r^D,D,---^{x). 



Including up to two derivatives, these operators give 
access to almost all continuum J^^ with J < 3. Suit- 
able linear combinations of these operators can be con- 
structed that transform as the irreducible representa- 
tions o f Table I. These are re lated to the operators 
used in iLiao and Mankel 120021 . 



Once a matrix of correlators, Cij(t)^ has been com- 
puted (for a given irrep), the mass spectrum follows 
from solution of a generalized eigenvalue problem that 
can be shown to be the quantum mechanical varia- 
tional solution. We solve 



C(t)v^ = K{t)C{tQ)v^, 



(2) 



for the eigenvalues Xaif) which are related to state 
masses, and for the eigenvectors v^ which are related 
to the overlap of our operators onto the mass eigen- 
states, the Zf in eqn ([l]). 

We computed correlators on quenched anisotropic 
lattices with a.s ^ 0.1 fm and o^^ ^ 6 GeV. Full de- 



2007 



tails can be found in |Dudek et al.| ^ ^ 

We show in figure [IJ the mass spectrum extracted 
for negative parity and charge conjugation. In the 
Ti representation we see precisely the level structure 
we expected, namely a ground state and three closely 
spaced excited states above. Looking at the states 
in the other irreps we see that one possible contin- 
uum spin assignment of the states in the first excited 
"band" would be to have two spin-1 states, one spin-2 
state and one spin-3 state^. We gain a good deal of 



■"^In the continuum the appearance of a e.g. spin-2 state 



Figure 1: Extracted mass spectrum for PC — listed 

by latti ce irreducib le representation. Operator labels 
listed in pudek et"ai.| [2007] 
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Figure 2: Extracted mass spectrum for PC = listed 

by assigned continuum spin. 



support for this hypothesis from studying the eigen- 
vectors extracted from eqn ([2|. Consider the lattice 
irrep projections of the "ai x V" operator: 



O^ 



ijk 



(9^. = o^-^'^ 






where |e*-^^|, Q*-^^ are Clebsch-Gordan coefficients for 
the lattice cubic group. In the continuum we know 
the form that the overlap of these operators onto a 
spin-2 state takes, so that 



iO^|2- 



■(p,r)) 



,ijk 



^jk (P^r), 



in both the three-dimensional T2 and the two-dimensional E 
corresponds to the five spin projections of a spin-2 meson 
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where Z is common to both. Z can be extracted from 
the eigenvectors and if it is found to be close in value 
in the T2 and E cases then we conclude that it is likely 
that we have a spin- 2 state. 

We apply this eigenvector inspection method wher- 
ever possible and where the result is conclusive we 
assign the continuum spin shown by the color cod- 
ing in figure [l] For the states above the first excited 
band this method gave inconclusive answers and for 
this reason we do not try to assign a continuum spin. 

This first calculation was performed only at one 
(quenched) lattice spacing and consequently our re- 
sults are not extrapolated to the continuum. Never- 
theless we present our results for continuum spin as- 
signed states in figure |2] along with experimental state 
masses taken from the PDG(Yaoet al. [2006 ) and po- 
tential model masses taken from iBarnes et al., [2005^ . 

It is clear that we are in agreement with the gross 
structure predicted by potential models, and in par- 
ticular we appear to have successfully extracted some- 
thing like the '^ (3686) /^^ (3770) system. We believe 
that this has not been achieved before in a lattice 
calculation. Extracted state masses appear to be sys- 
tematically high with respect to potential models and 
experiment - our suspicion is that this is due to some 
combination of computation at finite lattice spacing 
and the quenched approximation^ - this hypothesis 
can be tested with further calculation now that this 
method has been demonstrated. 

Other PC combinations were also considered. In 
figure [3] we show our results for J++. It is clear that 
again we are observing masses systematically higher 
than the potential model states. That we miss the 
spin-4 state near 4 GeV may be related to the fact that 
our operators, which have a maximum of two spatial 
derivatives do not have any overlap with spin-4 mesons 
in the continuum limit. This could be remedied by 
enlarging the operator basis. 

With PC = H — , as well as spin-singlets with odd 
J, one also has the possibility of exotic quantum num- 
bers, i.e. those not accessible to a qq Fock state. In a 
quenched heavy-quark calculation these can only arise 
through non-trivial gluonic excitation giving rise to 
states usually described as "hybrids". Our extracted 
mass spectrum is shown in figure |4] where exotic states 
with 0+~, 2^ quantum numbers appear above 4.5 
GeV. 

With PC — — h, the odd-J states are exotic. Our 
extracted mass spectrum listed by lattice irrep is 
shown in figure [5] This case demonstrates the dif- 
ficulty in continuum spin assignment; the set of five 



px V t 



pxV^ 
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J++ 



^In particular the problem of scale setting when one has an 
incorrect running of the coupling. In |Dudek et al.| [2007| the 
effect of finite box size was tested and was found not to be the 
source of the level raising effect. 
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Figure 3: Extracted mass spectrum for PC = -|--|- listed 
by lattice irreducible representation and continuum spin 
assigned states. 



levels near 4.3 GeV could, on the basis of their mass 
degeneracy, be interpreted either as a single ^ and 
a single (non-exotic) 4 ^ or as two ^ states, an 
exotic 1 ^ and a 2 ^ . In previous cases we used the 
eigenvector inspection method to break these ambi- 
guities, but unfortunately here the method produces 
inconclusive results. We display the two possible spec- 
tra in figure [6] where we note that the potential model 
does have a 4 ^ state in this mass range. 

It is worth pointing out that previous studies of 
the 1 ^ state in charmonia have not taken into ac- 
count the spin ambiguity and hence they may have 
in fact reported the mass of a non-exotic 4 ^ state. 
It is clear that further study with more operators and 
higher statistics is needed in order to make a definitive 
statement. 

We believe that we have demonstrated the power 
of using a variational solution in a large, carefully 
constructed operator basis to extract excited states 
in lattice QCD. Of course there remain numerous is- 
sues to deal with, including the effect of multiparticle 
(DD) states when one relaxes the quenched approxi- 
mation, but given that they too are orthogonal states 
we should be well-equipped with the method outlined 



74 8 



Proceedings of the CHARM 2007 Workshop, Ithaca, NY, August 5-8, 2007 



TT X V sm Q0(2) X D s: 

bl X D sm bi X P an 

ai X IB gjf, tti X IB ^!: 



bl X D s: 

^1 ^ J^ ^''' 



J+- 



:^^ 














1 ^=0 1 




1 J=^ 1 




1 


I 


1 






^ ^^^4515±46 : 










1 J = 2 1 






4553±14 




1 ^-3 1 


1 ^^ r 




^ 








1 1 


4166±70 


J=? 






3523±6 





1, 3, 4 ... 2, 3, 4 ... 





- 














5000 


\expt(PDG)] 


1 lanice i 




e: 


1 
















4500 


_ ^r~~ 




ICZ 


1 










- 


451 


±56 




1 


1 












4127±39 


1 q.m. P-wave \ 












4139±42 


j^^* 


^^^^^ 




4000 
















3000 









0+-(/ico) l+-(/ic) 2+-(/ic2) 3+-(/ic3) 4+-(hc4) 

qq exotic qq exotic qq exotic 



Figure 4: Extracted mass spectrum for PC = H — listed 
by lattice irreducible representation and continuum spin 
assigned states. 



75 am T47S sm *^0(2) X V sin TT X D s, 

^1 X V ,"'i bi X V am bi X V a 

p X B ,^^, p X B ^^, p X B s, 



bl X V a: 

p X B a",: 



- 


r 


4695±105 


4926±51 






1 J-o 1 




1 ^-' 1 










1 J = 2 1 










4379±23 


4417±60 


1 J-3 1 


- 


4330±61 


4327±36 


J=? 


1 


1 1 








1 






3876±44 










- 


3691 ±24 




" 


3023±2 





Figure 5: Extracted mass spectrum for PC = — h listed 
by lattice irreducible representation. 
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Figure 6: Two possible continuum spin interpretations of 
extracted mass spectrum for PC = — h- 



Proceedings of the CHARM 2007 Workshop, Ithaca, NY, August 5-8, 2007 



3. Radiative transitions 



Sub-I)5-threshold charmonia have very narrow 
widths such that radiative transitions between them 
constitute considerable branching fractions, these 
have been measured by a range of experiments and 
their relative magnitudes give us clues to internal 
structure. 

These transition widths can be computed using lat- 
tice QCD by considering three-point correlators of the 
type 

C^^J(t) ^ {O^{x,tf)V^{y,t)Oj{6,0)), 

where Oij are operators having overlap with meson 
states and V^ is a lattice representation of the vector 
current. For example we could extract the J/ip -^ r]cl 
matrix element from the large Euclidean times value 
of the correlator 



X 

X ( [^75V] {x, tf) [V^T^^] (y, t) [V^T.V'] (0, 0) 



Correlators of this type, using only point-like opera- 
tors were e valuated (details can be found in Dudek 
et al.|[2QQ6] ) and transition form- factors extracted for 
a set of transitions between J^^ ground states. 

The J /il) -^ rjcj transition form-factor shown in fig- 
ure [7| is the most statistically precise signal, but it suf- 
fers from a large systematic issue related to quenching. 
It is well known that the experimental hyperfine split- 
ting in charmonium is not reproduced well by stud- 
ies utilizing the quenched approximation. As such we 
have an ambiguity when computing the phase space 
that is required to scale a matrix-element to a width 
(or vice- versa) - should we use the experimental value 
or the value extracted from the spectrum portion of 
our lattice calculation? In figure [7| we show the ex- 
perimental width^ scaled to a matrix-element by both 
possibilities and the lattice data fitted with an expo- 
nential in photon virtuality, Q^, used to extrapolate 
back to Q2 = 0. 

A transition with reasonable statistical precision 
and a very small phase-space ambiguity is the electric 
dipole transition XcO -^ J/i^l- Our results are shown 
in figure |8] where the fit uses a form motivated by the 
quark model. Note the points at slightly timelike Q^ 
are not included in the fit - the agreement with the 
extrapolated curve then lends support to the fitting 
form used. 

Results fo r other transitions can be found in lDudekl 
et al. |2006| as can comparison of the lattice results to 
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Figure 7: Transition form factor for J/ip -^ rjc'^. 
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Figure 8: Electric dipole transition form factor for 
XcO -^ J/^1- Experimen tal data from PDG (2005) and 
CLEO dAda^ et al.| [2005] ) 



^We note that there is ongoing work at CLEO to confirm 
the single measurement from Crystal Ball 



potential model expectations. Work is currently un- 
derway combining the excited state technology of the 
first section with the radiative transition technology to 
make it possible to study transitions involving excited 
and high-spin states. This would include experimen- 
tally measured transitions like '0(3686) -^ XcJl- 
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4. Two-photon decays 



At first sight it is not clear how one would go about 
evaluating the matrix element for the process rjc -^ 77 
in lattice QCD. In the previous section we outlined 
how to extract the matrix element for a radiative 
transition between two QCD eigenstates from a three- 
point function evaluated at large Euclidean times. 
This issue here is that the photon is not an eigen- 
state of QCD - taking a vector interpolating field to 
large Euclidean time would not yield a photon state, 
but instead the lightest QCD vector eigenstate (the 
J /ijj in this case). 

However, all is not lost, for while the photon is not 
a QCD eigenstate, it can be constructed from a linear 
superposition of QCD eigenstates. The precise field- 
theoretic mechanism for this is the LSZ reduction. 
The connection in Euclidean space-time, for a differ- 



ent physical process, is made in Ji and Jung 2001 



and for the process in question an outline appears in 
Dudek and Edwards [2006 . The end result is that 
the following relationship connects the matrix element 
of interest to a Euclidean three-point function com- 
putable on the lattice: {ric{p)\^{qi^\i)^{q2-,^2)) ^ 



x( / d^xe-'P-^0{x,tf) / d^ye''i'^-yV''{y,t)V''{0,ti) 



(3) 



The difference with respect to the radiative tran- 
sitions between hadrons considered above is that an 
integral over the Euclidean time position of a vector 
source is now involved. 

The details of the lattice computation of th is object 
can be found in iDudek and Edwardsl [20061 , here we 



mention only that an isotropic lattice was used. In fig- 
ure^a) we display the integrand of equationpj having 
computed with an operator '^75'^ fixed at tf = 37, a 
conserved vector current insertion at t = 4, 16, 32 and 
a vector interpolating field at all possible source posi- 
tions, ti = ^ 37. It is clear that provided one is not 
too close to the dirichlet wall or to the sink position, 
one can capture the entire integral by summing times- 
lices. In figure loFb) the results of summing timeslices 
to compute the integral for all possible insertion posi- 
tions and a number of Q^ are shown - clear plateaus 
are visible at intermediate times indicating dominance 
of the r]c over the possible excited states. 

Given the confidence that the integral can be cap- 
tured on a lattice of this temporal length, one can 
use a much faster method to compute the transition 
form-factor that places the sum over timeslices into 
a "sequential source" , reducing the computation time 
by a factor of 0{Lt). Results using this method are 
shown in figu re p^O] along with PD G values and results 
inferred fromlUehara et al. 2007 . 
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Figure 9: (a) Integrand in equation [3] at three values of 
vector current insertion time (t = 4, 16, 32) with pseu- 
doscalar sequential source at sink position tf = 37. (b) 
Pseudoscalar two-photon form-factor as a function of time 
slice, t, from equation [3] First six time slices ghosted out 
due to the Dirichlet wall truncating the integral. 



Of course here the errors displayed on the lattice 
data are statistical only and must be augmented by 
an uncertainty due to scaling from our fixed lattice 
spacing to the continuum and one related to the lack 
of light-quark loops within the quenched approxima- 
tion. This is the first demonstration of this method, 
such controlled studies will doubtless follow now that 
efficacy has been demonstrated. 



5. Summary 

Several new techniques have emerged that much ex- 
pand the range of charmonium quantities that can 
be considered in lattice QCD. Initial studies with 
quenched lattices are clearly systematics dominated, 
but this can be expected to be improved in the near 
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Figure 10: (a) r]c -^ 77* amplitude, (b) Xco -^ 77* ampli- 
tude. Fits are one-pole forms as described in |Dudek and| 
Edwards! [20061. 



future by use of dynamical lattices, in particular the 
anisotropic dynamical lattices being generated under 
USQCD at Jefferson Lab. These same methods ap- 
plied to the light quark sector will provide invaluable 
information for future meson spectroscopy projects 
like GlueX. 
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